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The concept of "broken symmetry", that the 
symmetry of the vacuum may be lower than the 
Hamiltonian of a quantum theory, plays an im- 
portant role in modern physics. A manifestation 
of this phenomena is the Higgs boson in parti- 
cle physics 1 whose long awaited discovery is im- 
minent. An equivalent mode in superconductors 
is implicit in the early theories of their collec- 
tive fluctuations^. Spurred by some mysterious 
experimental results, 4 the theory of the oscilla- 
tion of the amplitude of superconductivity order 
parameter, which is the equivalent to the Higgs 
modes in s-wave superconductors 5 and its iden- 
tification in the experiments, was explicitly pro- 
vided^. It was also shown that a necessary con- 
dition for this to occur 7 is the emergent Lorentz 
invariance in the superconducting state while the 
metallic state and the region just below T c is man- 
ifestly non-Lorentz invariant. Here we show that 
d-wave superconductors, such as the high tem- 
perature Cuprate superconductors, should have a 
rich assortment of Higgs bosons, each in a differ- 
ent irreducible representation of the point-group 
symmetries of the lattice. We also show that 
these modes have a characteristic singular spec- 
tral structure which can be discovered in Raman 
scattering experiments. 

The order parameter \I/ in s-wave superfluids and su- 
perconductors is a complex number. The oscillation of 
the phase of ^ is the Bogolubov (Nambu-Goldstone) 
mode, which is massless at long wavelengths in a neutral 
superfluid. This phase mode may be understood as the 
azimuthal oscillation of a particle near the bottom of a 
Mexican hat potential, depicted in Fig. In a charged 
superconductor, it moves to the frequency of the plasmon 
in a gauge invariant theory coupling phase modes to elec- 
tromagnetism 8 . While very interesting as the VK-boson 
in particle physics, it tells nothing new about excitations 
of the superconducting state. In contrast, the amplitude 
mode, which oscillates in the radial direction does not 
couple to charge and has an excitation gap (mass gap) 
at long wavelengths equal to twice the superconducting 
gap 2 A. This is just where the continuum of particle-hole 
excitations begin, hence it is heavily damped and usually 
unobservable. Special situations which lower its energy 
are therefore required to detect this mode in s-wave su- 
perconductors P Interesting related modes have also been 
discussed in superfluid 3 H&^. 

In this Letter we show that superconductors with lower 
symmetries support additional amplitude or Higgs modes 
labeled by the point group symmetry in which the defor- 
mation of the order parameter occurs. As expected, one 




FIG. 1. Pictorial representation of the effective potential cor- 
responding for a) s-wave Higgs mode and b) additional non s- 
wave Higgs modes. The Nambu-Goldsone mode (black circle) 
oscillates in the azimuthal direction, whereas the amplitude 
(Higgs) mode (red curve) oscillates in the radial direction. In 
(a) the constant curvature of the effective potential results 
in an angle-independent finite mass gap. In (b) the effec- 
tive potential plotted for the non-conventional Higgs modes 
V(pi,6i) = (a + b sin 2 #)(p 2 — A 2 ) 2 exhibits two- fold symme- 
try leading to a periodic angular dependent curvature. The 
non s-wave amplitude mode results in a more massive fluctu- 
ation at Oi = 7r/2 than Oi — 0, in contrast with the s-wave 
Higgs mode. This leads to an angular dependent energy uj{0i) 
with the minimum and maximum occuring at Oi — 0, 2tt and 
Oi — 7r/2, 3tv/2 respectively. 



of these modes is the conventional s-wave Higgs mode 
which appears at 2A and is likely to be overdamped. 6 
The other amplitude modes, however in general, have 
lower energies than 2A, where the damping is smaller. 
They may therefore be more easily identified. As these 
additional modes correspond to deformation of the or- 
dered state to different irreducible representations with 
variable relative phases Oi, the Higgs modes acquire a 
characteristic singular two-peak lineshape which is de- 
rived here. 

Besides the U(l) gauge symmetry, anisotropic su- 
perconductors are also invariant under a point group 
symmetry determined by the crystal lattice structured 
For definiteness, we consider a two-dimensional uncon- 
ventional d-wave superconductor with d x i_ y i ordering on 
a square lattice, point group symmetry The high 

temperature Cuprate superconductors belong to this cat- 
egory. The nature of the additional amplitude modes is 



sketched in Fig. [2] and corresponds to excited states with 
admixtures of additional d x 2_ y 2, d^-wave, g xy (x 2 -y 2 )- 
wave and s + 2 -wave components to the ground 

state, respectively. 

We represent the ground state and the oscillations 
about it by the order parameter, 



tt(Q,k) = tf (k)+<S*(Q,k,w)e 



iO(Q) 



(1) 



^o(k) is the uniform ground state which we assume 
to be in the Bi g , i.e (k 2 — k 2 ) symmetry with phase 
= 0, £\I/(Q,k, uo) are the amplitude of the devia- 
tions representing the collective modes with total cen- 
ter of mass momentum Q and internal momentum k 
with phase #(Q). At long wavelengths, 5^ may be writ- 
ten as a separable function of Q and k. The k de- 
pendence is expressed in the four one dimensional even 
parity irreducible representations (B\ gi A\ gi £?2#, ^2g) of 
the D/t point group symmetry. For notational simplic- 
ity, we will represent £\I/(0,k) as linear combinations of 
(j)i(k) = |^i(k)| exp (iOi); i = 0, 1, 2, 3, respectively. In the 
limit Q = 0, the field theory is given by the Lagrangian 
(see supplemental material), 



C = ^\d t &\ 2 +a z \&\ 2 -b z \&\ 4 



(2) 
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i<j 

We include only second order time- derivatives; this is 
only valid well below the Ginzburg-Landau regime near 
T c , where a first derivative in time representing dissipa- 
tion dominates and Higgs mode cannot occur due to lack 
of Lorentz invarianceA We have introduced two distinct 
set of parameters Cij and dij in ([2| so that the energy of 
the collective modes depends on the relative phase 0i be- 
tween the assumed ground state representation and the 
others. This is required by symmetry and introduces dis- 
tinctive features in the spectra of the collective modes as 
we see below. 

The equations of motion using ([2| give the energy of 
the collective modes at Q = to be 



WiWi) = ±V feo + *o sin^))|^o| 2 + a,, (3) 

here 0i is the relative phase of the i / order param- 
eters with respect to the ground state order parameter 
|^o | 2 = — ao/2&o- The (k 2 — k 2 ) order parameter as- 
sumed for implies ao < for T < T§ = T c (where 
T c is the critical temperature). a[s(i ^ 0) remain pos- 
itive as T approaches T c from below, > d^ > 0's 
are expected because of the co mpetition between differ- 
ent order parameters, ujq = y4^oT^oP corresponds to 
the simple s-wave Higgs mode of the d-wave supercon- 
ductor and appears at 2 A. The energies uoi correspond 
to fluctuations of the d k 2_ k 2 order parameter in which it 
deforms to other point group symmetries as depicted in 
Fig. |2] 
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FIG. 2. Pictorial representation of the additional Higgs or 
amplitude modes of the d-wave superconducting order pa- 
rameter predicted in this Letter. Each mode can be labeled 
by an irreducible representation of the point-group symmetry 
of the lattice in which the deformation of the order parame- 
ter occurs (see text for details). For the case of d k 2_ k 2 order 
parameter depicted above these amplitude fluctuations are 
different admixtures of d x 2_ y 2-w&ve ("breathing mode"), d xy - 
wave ("rotating mode"), g xy ( x 2_ y 2ywa,ve ("clapping mode") 
and s + ^( x 2_ ;y 2)2-wave ("osculating mode") components to 
the ground state, labeled from top to bottom. 



The mass uji of the modes can be estimated from 
general considerations and by comparison with s-wave 
Higgs mode. In order to compare the energies uoi 
with ujq one can gain insight by using a two-parameter 
Landau-Ginzburg energy functional in the parameter 
subspace (<fo,^i). The phase diagram in this subspace 
allows for three broken symmetry phases a) |</>o| 2 = 
-a /(2&o), |^i| 2 = 0, for a < 0, a* > ; b) |0 O | 2 = 
0; \(j)i\ 2 = —a,i/(2bi) for ao > 0, ai < 0; and a mixed phase 
c) |(/>o| 2 7^ 0, \(j)i\ 2 ^ which only appears for a^ < and 
ao < 0. Since we assume that the broken symmetry su- 
perconducting state has d k 2_ k 2 order, we must require 
that |ao| > \di\ for a^ < and ao < 0. In order to avoid 
a second order transition to the mixed phase we must 
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satisfy 



Cio<2y/bobi, and c i0 > 2b - 



(4) 



which establishes an upper and a lower bound on the 
energies coj. 

In order to estimate the values for bi and we assume 
an attractive potential is dominant for all the irreducible 
representations, 

V(k - k') =V! + V a (k)a (k') + V^a 2 {k)a 2 {k r ) 

+ V4ao(k)a 2 (k)ao(k , )a2(k / ), (5) 



k 2 and a 2 (k) 



k x k y , with Vi < 



where a (k) ,. 

for all values of i and |Vb| > |T^|(z + 0). This is 
a natural assumption for the d xy symmetry, since the 
difference from d x 2_ y 2 arises only due to the anisotropy 
in the density of states, and so also for the d xy [ x 2_ y 2^ 
case. No such strong argument can be given for the s- 
wave case and so a repulsive potential is allowed for this 
caseP^In the "weak coupling" limit NpVo << 1, where 
Np is the density of states evaluated at the Fermi energy, 
an estimate for the values of and bi gives 



a 



Vi \og{T/T?) 
V (T-T c ) 



bo 



(6) 



for T ~ T c . Since a { > for T - T c the energies ^ 
starting initially at a non-zero value decrease in magni- 
tude for temperatures below the transition temperature 
T c as di — >• for T — >> Tf, whereas co>o increases in mag- 
nitude as the temperature is lowered from the transition 
temperature T c . This can be seen from a combination of 
energy expression ([3| and the upper-bound on the val- 
ues of Cio ~ dio which gives Qo < 2bo(Vi/Vo) 1 ^ 2 << 46o- 
This indicates that there exist temperatures T* < T c 
where uji « cjo, thus establishing an upper bound on 
the energies c^. 

The lower bound in Eq. [4] follows from the condition 
that energies uji are always positive, so that no transition 
from the chosen ordered phase is allowed. The collective 
mode energies uoi for the Higgs modes as a function of 
a combination of the phenomenological parameters and 
superfluid density are depicted schematically in Fig. [3^i. 
An examination of Eq. ([3| and subsequent considera- 
tions reveals that uoi at T = for (i ^ 0) is simply the 
difference of the ground state energy of the i-the symme- 
try from that of the realized (i = 0) symmetry, as could 
have been guessed at the outset. 

Since, at low energies, the quasi-particle density of 
states in a d-wave superconductor is proportional to the 
energy, the lower the energy of the modes in Fig. [3^i, the 
less they are damped. All the Higgs modes in d-wave su- 
perconductors, being oscillations of the amplitude of the 
superconducting condensate are neutral spin modes. As 
such they do not couple to the usual external probes. In 
the case of s-wave superconductors^, they could be dis- 
covered only through appearing in the self-energy of the 
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FIG. 3. (a) Schematic behavior of the energies for the ad- 
ditional Higgs modes at ft = as a function of the cou- 
pling constants and superfluid density cn/cio\^o\ 2 at T = 0. 
The modes are labeled by their Raman scattering geometry, 
the black (solid) line corresponds to the rotationally symmet- 
ric Aip-Higgs mode, whereas the red(dotted), blue(dashed) 
and green (dot-dashed) lines correspond to the additional Ai 9 - 
,B2 9 -,Bi g - Higgs modes(see text for details). We have assumed 
an attractive potential in all the irreducible representations, 
(b) Schematic representation of the line shape associated to 
the "clapping" i^-Higgs mode (chosen arbitrarily) indicat- 
ing the energy continuum and square-root singularities at the 
edges of the energy spectrum 003(0) and cj3(vr/2). The con- 
ventional Ai^-Higgs mode (which is likely to be overdamped) 
is shown as a reference. 



superconducting state of phonons which promotes super- 
conductivity, and steal intensity from them. Similarly, 
for cuprates, we expect that if the broad quantum-critical 
fluctuations, whose q — >• limit is visible in Raman scat- 
tering, promote super conduct ivit}B^, will partially give 
its weight to the Higgs modes. Elementary considera- 
tions indicate that uj 2 or the breathing mode (which is 
likely to be over damped) occurs in the s-wave or A\ g 
symmetry because for it 5^/ also has (x 2 — y 2 ) symmetry, 
the rotating mode occurs in the A 2g symmetry because 
for it 5^ has xy symmetry, the clapping mode occurs 
in the B 2g symmetry because for it Sty has xy(x 2 — y 2 ) 
symmetry and the osculating mode occurs in the B\ g 
symmetry because for it Sty has s-wave symmetry. The 
line shapes, which can be calculated from Eq.[3J exhibit a 
two peak structure with square root singularities at the 
edges of the energy spectrum as shown in Fig. |3Jd (see 
supplementary material). The actual observation of the 
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Ai g mode may occur as a sharp peak below 2A through 
coupling to the continuum. Indeed a mysterious intense 
mode in the A\ g channel has already been detected 14 . 

In the supplemental material we deduce the La- 
grangian for the gradient terms to derive the leading Q 
dependence of the energies. We find that there are inter- 
esting couplings between phase and amplitude modest 
which are not influenced by Coulomb interactions, be- 
cause of counterflow in the excited states of supercur- 
rents in two different symmetries keep the system charge 
neutral. However, to quadratic order in Q, the energy 
spectrum oJi(Oi) remains unchanged due to this coupling, 
only acquiring a quadratic dependence in the wavevector 
Q (see Eq. (10) in the supplementary material). Any 
effects of this coupling appear beyond quadratic order in 

Q 

The low-energy physics of many condensed matter sys- 



tems (lattice bosons near a Mott transition 16 , antiferro- 
magnet^, incommensurate charge-density wave and su- 
perconductors 3 ), close to a quantum critical point, is cap- 
tured by a Lorentz invariant critical theory 1 ^. One con- 
sequence of spontaneous breaking of a continuous sym- 
metry in a Lorentz invariant theory is the appearance of 
amplitude fluctuation or "Higgs" modes. We show that 
when this symmetry is endowed with an additional dis- 
crete space symmetry, a rich assortment of Higgs modes 
should be present. The number of these Higgs modes 
should be equal to the number of irreducible represen- 
tations of the discrete point group symmetry consistent 
with any internal symmetries (such as spin or valley) of 
the system. 
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